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Abstract. We prove that given a stress-free elastic body tliere exists, for sufficiently 
small values of the gravitational constant, a unique static solution of the Einstein equations 
coupled to the equations of relativistic elasticity. The solution constructed is a small defor- 
mation of the relaxed configuration. This result yields the first proof of existence of static 
solutions of the Einstein equations without symmetries. 



1. Introduction 

General relativistic effects generated by compact, isolated bodies, such as stars and 
even satellites, are of increasing importance in observational astronomy and experimental 
general relativity. Considering this fact, it is remarkable how little is known about solutions 
of the Einstein field equations for systems with spatially compact sources. The situation is 
not much better if we describe gravity by Newton's theory. 

The present paper provides the first existence result for compact, isolated, static elastic 
bodies in Einstein's theory of gravity. With the notable exception of collisionless mat- 
ter, essentially all previous results concerning compact, isolated, self-gravitating bodies 
deal with static or stationary fluid bodies. Under reasonable conditions, static fluid bodies 
are spherically symmetric, while stationary fluid bodies are axi-symmetric. Although the 
variational formulation of elasticity has strong similarities with that of fluid models, static 
elastic bodies may, in contrast to static fluid bodies, be non-symmetric. In fact, in this 
paper we prove, for the first time, existence of static solutions of the Einstein equations 
without symmetries. 

1.1. Compact bodies. Fluids and dust (i.e. a pressure-less fluid) are the conceptually 
simplest and most commonly used matter models. For a self-gravitating compact body, it 
is necessary to consider a free boundary problem with zero traction on the boundary. It 
was only recently that an existence proof was given by Lindblad for the Cauchy problem 
for a nonrelativistic perfect fluid with free boundary, in the absence of gravity |24, 23 1. For 
self-gravitating fluids, no results of this generality are known, and it is only in the static or 
stationary cases that results are available. 

If we assume that spacetime is static, the standard conjecture is that any isolated self- 
gravitating body, consisting of a perfect fluid, is spherically symmetric. This is known to be 
the case in Newton's theory, for a general equation of state. In Einstein's theory there is the 
work by Beig and Simon |8| which solves the problem for a large class equation of states; 
a proof in the general case remains to be found. For the case of stationary spacetimes in 
Einstein gravity, i.e. spacetimes with a timelike, non-hypersurface orthogonal Killing field, 
it is known under certain additional assumptions on the thermodynamic properties of the 
fluid that there exists an additional, rotational Killing field, so that stationary spacetimes 
containing fluid bodies are axisymmetric ||251 . 
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In spite of the symmetry restrictions discussed above, there are rich classes of station- 
ary and static solutions describing isolated bodies, even in Newtonian gravity. The almost 
completely forgotten work of Leon Lichtenstein from roughly the period 1910-1933, pro- 
vides existence results in the Newtonian case for rotating fluid solutions in various con- 
figurations, see the book ll22ll . Inspired by Lichtenstein, Uwe HeiUg showed in 1995 the 
existence of stationary rotating fluid solutions in Einstein's theory (15). 

These results allow one to construct stationary fluid solutions with slowly rotating al- 
most spherical balls. Further, one has rings, rings around balls and families of nested rings. 
It seems almost impossible to get an overview on all possibilities. New solutions can often 
be constructed as perturbations of known solutions. For example, starting with a static, 
spherically symmetric fluid ball whose existence can be shown by using ODE techiques, 
one may prove the existence of a rotating solution with small angular velocity. This is 
essentially what was done by Lichtenstein and Heilig. 

The Vlasov matter model is a statistical description of weakly interacting particles. It is 
conceptually more difficult to work with than fluids but has been used very successfully in 
various circumstances. For a survey of known results, see Ii28il. The existence of various 
dynamical and time independent solution has been demonstrated. All the known stationary 
and static solutions have axial symmetry. 

1 .2. Elasticity and relativity. Elasticity is of course one of the oldest topics of theoretical 
physics, with origins that can be traced back to the 17'th century. The book by Marsden 
and Hughes |26| gives a modern treatment of elasticity. Already in 1911, Herglotz 1161 
gave a formulation of elasticity in special relativity. There are various formulations of 
elasticity in the framework of general relativity, see for example Rayner 1271 . Carter and 
Quintana 111], Kijowski and Magli ifTSl [191 , Christodoulou 1121 [131 . to name just a few 
important works. Strangely enough the problem of existence of static or dynamical self- 
gravitating elastic bodies in Einstein's or Newton's theory of gravity has, to the best of our 
knowledge, until recently not been considered. The only exception is for the spherically 
symmetric case. Even in non-relativistic elasticity quite little is known. The first existence 
theorem in three dimensional static elasticity was given by Stoppelli in 1954 f29^. 

About six years ago, two of the present authors (R.B. and B.S.), motivated by this state 
of affairs, initiated a program to develop existence results for elasticity in the setting of 
Einstein's theory of gravity. We first showed the existence of static solutions describing 
elastic bodies deformed by their own Newtonian gravitational field |5|, and later estab- 
lished the existence of relativistic elastic bodies deformed under rigid rotation jSj. For 
time dependent solutions local in time and space (no boundary conditions) uniqueness was 
proved in |[T3l . and for existence see H. In 19] an existence theorem for the motion of a 
free elastic body in special relativity is given. 

Elasticity can be described as a Lagrangian field theory l4l [T3]| . and hence the action 
for self-gravitating elastic bodies is derived by simply adding the gravitational Lagrangian 
to the Lagrangian for elasticity. The basic matter field is a map, the configuration, from a 
region in spacetime to the body, an abstract 3-manifold whose points label the constituents 
of the elastic body moving in spacetime. The stress of a configuration is determined by 
the stored energy function which completely fixes the matter model. This formulation 
can be used in a non relativistic spacetime, in special relativity and in Einstein's theory. 
Elastic materials where the stress is determined by a stored energy function are usually 
called hyperelastic. Diffeomorphism invariance, a necessary condition in Einstein's theory, 
implies that the stored energy function satisfies the additional condition known in the non- 
relativistic case as material frame indifference. 
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Once the stored energy function is given, the variational problem as well as the Euler- 
Lagrange equations are determined. In particular, one is led to consider Einstein's field 
equations with an energy momentum tensor which is determined by the deformation. The 
elasticity equations are a consequence of the conservation law in Einstein's theory. 

As the configuration is a map from spacetime to the body, we have a free boundary value 
problem. To deal with this difficulty, one reformulates elasticity using deformations, i.e. 
maps from body to spacetime, as the basic variable. In this setting, known as the material 
picture, one has a fixed apriori known boundary. 

Let us now consider the static self-gravitating bodies in general relativity. In this case, 
the theory can be given a a variational formulation on the quotient space of the timelike 
Killing field. Thus, in order to construct a static self-gravitating body in Einstein gravity we 
start from a relaxed elastic body without gravitational field and determine the deformation 
of such a body under its own gravitational field for small values of the gravitational constant 
G. To do this, it is convenient to choose a stored energy function for which there exists a 
configuration which is stress free, i.e. which satisfies, together with the Minkowski metric, 
the Einstein field equations for G = 0, as well as the the elasticity equation. We start from 
this background solution and construct nearby self-gravitating solutions for small values 
of G. 

As mentioned above, self-gravitating static fluid bodies in general relativity are known 
to be spherically symmetric for a large class of equations of state. In fact, the result proved 
in this paper provides the first example of a static solution to Einstein's field equations 
which is not spherically symmetric. 

It is worth pointing out that the approach used in this paper cannot be applied to the 
problem of constructing for example an elastic neutron star since in this case, there is no 
nearby stress free configuration. To deal with this problem one would have to choose a 
stored energy function in which the shear strains are much smaller then the hydrostatic 
compression. Then one could begin with a spherically symmetric solution in which the 
radial pressure would be different from the tangential pressure. One could then use the 
methods developed in this paper to construct nearby solutions which are not spherically 
symmetric. 

Finally we remark that it may be argued that the result proved in this paper is very weak, 
since G is required to be sufficiently small. However, it should be noted that we make no 
restriction on the shape of the body. For example, one may consider two very large bodies 
connected by a very thin neck. In this situation it is clear that if we make gravity too strong 
(i.e. G too large), the neck wiU break and hence there can be no static solution for such a 
configuration for arbitrarily large values of G. Thus, without restrictions on the shape of 
the undeformed body we can not expect a stronger result. 

1.3. Overview of this paper. In section|2]we give some analytical preUminaries. We also 
for the convenience of the reader review some basic ideas from linearized elasticity which 
we will make use of. Further, we prove some results which will be used concerning Bianchi 
identities for weakly differentiable metrics, and concerning the divergence of tensors wich 
compact support. Section[3]presents the gravitational field equations in space and the elas- 
tic equations on the body. We use harmonic gauge to make the field equations elliptic. An 
important step is to extend the body to R'^, the extended body, and to extend the deforma- 
tion, the inverse configuration, to a map from the extended body to physical space, which 
in our case also is W". Then we move the field equations from space to the extended body. 
In this way we obtain a quasilinear system of partial differential equations, the reduced 
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Einstein equations in material form, were the geometrical unknowns are defined on the 
extended body and the elastic variables on the body. 

In section|4]we formulate the reduced Einstein equations in terms of a non-linear map- 
ping between Sobolev spaces and calculate the Frechet derivative of the map at the relaxed 
configuration. It contains esssentially linearized gravity and hnearized elasticity. The lin- 
earized operator is Fredholm with non trivial kernel and range. The geometric reason for 
this is the combination of diffeomorphism invariance of the Einstein equations and Eu- 
clidean invariance of the background solution. Using an approach to some extent inspired 
by II2TI . we define a projection operator such that we can use the implicit function theorem. 
This way we obtain for small G a solution of the reduced field equation together with the 
projected elasticity equations. 

Section |5l which is the heart of the paper, contains a proof that the solution to the 
reduced, projected system obtained using the implicit function theorem, is in fact a solution 
of the full system of equations for the self-gravitating elastic body. At first it might seem 
that there are two possibilities to prove this. On the one hand, if the exact elastic equations 
were satisfied, a standard argument using the Bianchi identities would imply that also the 
harmonicity condition is satisfied and we have in fact solved the full field equations. On 
the other hand, if we could show that the full Einstein equations hold, the Bianchi identities 
would imply we also have a solution to the exact elastic equations. 

In fact, none of these two alternatives are applicable, and one must prove both properties 
simultaneously. To do this a type of bootstrap argument must be used. It is worth men- 
tioning that the boundary condition of vanishing normal traction is essential. We could in 
principle solve the projected elasticity equation together with the reduced field equations 
for a boundary condition which prescribes the position of the boundary in space. However, 
it would then in general not be possible to show that all the Einstein field equations are 
satisfied. This is consistent with the fact that fixing the surface of a body in space is not a 
physical problem in Einstein's theory. 

The result of this paper is proved in a way which is completely different from the analo- 
gous result in our Newtonian paper O. In an appendix (AppendixlBli we add an outline of 
the proof of the Newtonian result, found and kindly communicated to us by an anonymous 
referee, which exactly follows the pattern of the present work. 

2. Preliminaries 

The following index conventions will be used. Upper case latin indices A, S, C, . . . 
take values 1,2,3, lower case latin indices i, j, fc, . . . take values 1,2,3, and greek indices 
a, /?, 7, . . . take values 0, 1, 2, 3. 

We will make use of Sobolev spaces W^'^ on domains and the trace spaces B^'P, as well 
as weighted Sobolev spaces W^'^. Unless there is room for confusion, the same notation 
will be used for spaces of tensors and vectors, as for spaces of scalar functions. The rest of 
this section collects some notions and facts from analysis which shall be needed. 

2.1. Sobolev spaces on domains. The books HKTO] are general references for the mate- 
rial discussed in this section. For an integer fc>0, l<p<oo and a domain Q,, W^''p{H) 
is the closure of C°° [il] in the norm 

\a\<k 

Further, we shall need the Nikol'skii-Besov spaces B^'P — B^p, which are the trace 
spaces for the Sobolev spaces. These are Banach spaces with norm defined on M" for 
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s > 0, 1 < p < oo, by 

||u||b=.. = II^IIl. + \h\-'--+^^)\\^lu\\l,dhj'' 

where a is the smallest integer strictly greater than s, and A/i is the difference operator. 
There are versions of many of the facts stated in this section also for p = I, and p = oo, 
see the references given above. Note B'^'P ^ w'^'P, except for the case p — 2. Let fc > 1. 
A well known fact is that for a domain ft C K" with C*^ boundary, then for 1 < p < oo, 
the trace tiga has the property 

We shall make use of this fact in the case fc = 1. Further, under these conditions, there is 
a bounded linear extension operator E : W^'P{^) -> W^'P{W) n C°°(R" \ Vt). Here 
= ri U dVL denotes the closure of Vl. In fact, for the last mentioned result to hold, it is 
sufficient to assume that Vl has a Lipschitz regular boundary. 

2.2. The boundary problem of linearized elasticity. The book f30l is a general refer- 
ence for the material discussed in this section. Let Ai^ k be a fourth order elasticity tensor 
on a domain Vl C M" , i.e. A has symmetries 

fc — k — -^i I — ^fc I 

Let 

cj{uy = A^k'diuK 

and defined the operator L by 

Lui = dja{u)i-' 

L is strongly elliptic if 

A^k'v'^jV'^i>0, foralU,r;e]R". 

For applications in elasticity it is natural to assume there is a positive constant A such that 
for all symmetric n x n matrices ip'^j 

(2.1) A|vp < < A-VP, 

see ll30l Chapter III] or ll26l S 4.3]. The pointwise stabiUty condition (12.1b implies strong 
ellipticity. We will assume that ( 12. Il l holds for the rest of this section. 
The Neumann type problem 

(2.2) Lui = hi, trdna{u)i^nj = n 
is equivalent to the statement that 



Let H'^{ri) = B'^'^{ri). Then A defines a bounded quadratic form on iJ^(fi). The 
radical Z of A is the space of ^ e iJi(fi) such that A{^,u) = for all u e H^{Q). 
It follows from (12.1b and the symmetry properties of A, that Z consists of all Euclidean 
Killing fields ^, i.e. fields of the form 

(2.4) C = a' + b'jx^ , a', constants, 6*^ = -Ir'i 



(2.3) = - / (f>% 

Jn 

for all (p S C°°{n), where A{u, v) is the symmetric bilinear form 

A{u,v) = / d.v'A.h'diu'' 
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Let {Q)e be the orthogonal complement of the radical, i.e. the space of u e ^ (ft) 
such that 

(2.5) / ^^u'^o, y^ez 

Under the above conditions, the quadratic form A{u, v) is coercive on H^{Vi)e. This fol- 
lows from the pointwise stability condition ( 12.11 ) and Korn's inequality, see ||30l p. 92]. 
Thus, by the Lax-Milgram theorem, we have that for any r^) g _ff (fi) x H~^/'^{d^) 
satisfying 

(2.6) / - / = 0, ve e Z, 

Jn Jdfi 

there is a unique u e which is a weak solution to ( 12.21 ). and which satisfies the 

estimate 

\\u\\h^ < C{\\b\\H-i(n) + l|T||ff-i/2(ao)) 

We will later refer to i2.6i as an equilibration condition. The physical meaning of the 
equilibration condition is that the total force and torque exerted by (6, r) is zero. Now, 
assuming dfl e (7*'+^, Aj-'fc' e u e W''+'^'P{n), one has from |2| an estimate of 

the form 

l|u||iyfe+2.p(n) < C{\\L{u)\\wk.p(n) + \\iidn<j{u) ■ n\\Bi-i/p.p(^on) + ||u||Lp(n)) 

Let M^'''P(0)e be the space of u e W'''P{n) such that the condition ( 12.51 ) holds. Then we 
have by the above that for b,T e W''^P{n) x B''+^-^^P^P{dn), satisfying ^M, there is 
a unique u € W''+'^^P{n)e which solves (O. In particulai-, in view of the above stated 
estimates, we have that the Hnear mapping W''+^-P{n) W''-P{n) x B''+^-^^P'P{dn) 
defined by 

(2.7) u {Lu, traoo-(u) • 12) 

is Fredholm. It follows from the discussion above that the cokernel of the operator defined 
by ( 12.71 ) is defined by ( 12.61 ). and that the kernel consists of Killing fields, of the form given 
in (1231). 



2.3. Weighted Sobolev spaces on M". The material which we shall need can be found in 

IS. Let n > 2, r = \x\, and let cr = (1 + r^)^/^. For k > 0, k integer, 1 < p < 00, S e R, 
define function spaces Wg''' as the closure of Cg° (M") in the norms 

We use the notation Lg for Wg'^. Decreasing S means faster decay. The spaces W^^^^^ are 
equivalent to the homogenous Sobolev space W'^'P without weight, in particular LP_^^^^ = 

LP. The dual space of Wg'^ is Wgi''^' with l/p+l/p' = 1, 6' ^ -6 - n. Weighted 
Holder spaces Cg'" can be defined analogously, see fJl. 

Weighted function spaces can be analyzed in terms of ordinary function spaces by the 
following standard construction. Let be a bump function with support in |x| < 2 and 
which equals 1 on |x| < 1. Define a cutoff function tp by ip{x) = 4){x) — (j){2x). Then 
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tp has support in the annulus 1/2 < |a;| < 2. For a function u, let uq = (pu, and let 
Ui{x) — ■ip{x)u{2'^x) be the dilatations of u. An equivalent norm for W^'P is given by 

oo 

Using this formulation, inequahties on compact domains may be systematically general- 
ized to weighted spaces. The following are some of the basic inequalities for the weighted 
spaces. 

(1) Inclusion: If pi < p2, S2 < Si, and u e L^^, then 

(2.8) My,^<C\\u\\^p 

(2) Sobolev I: \fn — kp> {) and p < q < np/ (n — kp), then if u G W^^'^, 

\\u\\l^^<C\\u\\^... 

(3) Sobolev II: If n < kp, then if m e Wg'^, 

\\u\\Lr <c\\u\\^.,., 

and in fact \u{x)\ = o{r^) as r cx). 

(4) Product estimate: If n < kp, and u G W^^^, v £ Wg^^, then with 6 = Si + S2, 

(2.9) ll"w||w'=.p < C'll"lln/'=.p|b||w'=.P 

Let A denote the Laplacian defined with respect to the Euclidean metric on M". We 
recall some facts about its mapping properties in the setting of weighted Sobolev spaces. 
Let E — {j : j integer , j 7^ 3 — n, . . . , — 1}. In particular, for n — 3, E consists of all 
integers. The elements of E are called exceptional weights. A weight 5 G M \ _E is called 
nonexceptional. Given (5 G M, define k~ (S) to be the largest exceptional weight < S. 

A basic fact is that A : Wg''' WgZ2"^ for fc > 2, 1 < p < 00, is Fredholm if and 
only if S is nonexceptional. In particular, for (5 G (2 — n, 0), A : Wg'^ WgZ-^'^ is an 
isomorphism. Let the operator L be of the form L ~ a^^ didj + h^di + c. Then with q > n, 
we will say that L is asymptotic to A, of order r < 0, if 

Note that the above conditions are stronger than those stated in |[3] Definition 1.5], and that 
we use the opposite sign convention for r. If L is asymptotic to A, then for 1 < p < q, 
L : Wg'^ L^g is bounded. 

The following version of elliptic regularity is easily proved using standard estimates for 
elliptic operators on domains (see for example lfT4l Chapters 8,9]) and scaling. Suppose 
L is asymptotic to A of order t < and suppose 2 < p < q. For u G Wg'^, such that 
Lu G i^_2' elliptic regularity gives u G Wg'^ and the inequality 

holds. See [3] Proposition 1.6] for a stronger version of elliptic regularity. 

If S is nonexceptional, there are constants C, R so that if m G Wg''^, the scale broken 
estimate 

\\u\\y^,2.p < C{\\Lu\\lp_^ + ||u||lp(b«)) 

holds, cf. Is] Theorem 1.10]. Here Bfl — {x : |a;| < R}. We will now state a consequence 
of this estimate which we shall make use of. Assume that S + t is non-exceptional. If 
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u G Lu G i^+T-_2^ then for exceptional values S + t < j < k (S), there are 

/jj (joo (gn^^ harmonic and homogenous of order j in M" \ Bfi, such that 

<5+r<j"<k-(i5) 

where v G Wg^^, and an estimate of the form 

holds. Here \ \h^ \j is a suitable norm of for homogenous, harmonic functions, for example 
\\h^\\j — ||r^^7i^||ioo(52^), with S2R = {x : \x\ = 2R}. In paiticular, if kerL = on 
Wg'^, then the above estimate takes the form 

To make the above explicit, suppose n — 3, —1 < S < and —2 < S + t < —1. Then 
k~((5) = — 1 and with u G Wg'^, Lu G -^5+^_2' ^^^^ 

u = — Q + V 
r 

where C is a cutoff function such that ( = 1 in R'^ \ _Bfl,and C = in and where 

V G Wg]^^ satisfies an estimate of the form 

\ci\ + \\v\\^2^.^<C{\\Lu\\i^.^^_^ + \\v\\bJ 

We shall make use of this estimate in section ISTTl 

For T < 0, p > n, define the space EI^'P of asymptotically Euclidean metrics on R" as 
the space of /i^ such that 

h,, - % G W^^P 

where 5ij denotes the flat Euclidean metric on R". Then EI^'P is a Banach manifold. 

Let Rij be the Ricci tensor of hij. We shall make use of the fact that if ft, G E^-p for 
p > n, the operators and 1-^ LVi = AhVi + RijV^ are asymptotic to A of order 
T. Here it should be noted that the principal part of the operator L is the scalar Laplacian, 
acting diagonally. The results concerning elliptic operators that we have stated in this 
subsection generalize immediately to the case of elliptic systems of diagonal form. 

2.4. Bianchi identity. The Bianchi identity for weakly regular Riemann spaces will play 
an important role in this paper, and therefore we give a proof of this fact below. The 
considerations in this subsection are local, and we work in local Sobolev spaces, denoted 
by Wilf. By definition, / G Wj^f if for each compact domain U C M, f G W''-p{U). 
Fori <p< cx), let g be the dual exponent of p, such that 1 /p + l/g = 1. For non-negative 
integers k, we define W^^^''' as the dual space to Wl^f. Then d°'f G W^^^'^ if / G Lf^^, 
where k — \a\. 

We make note of the following product estimates. Suppose that p > n. If u G W^f, v G 
W^^\ then uv G W^f. To see this, differentiate the product and use Sobolev imbedding. 
Further, if u G W^f, w G Wi^^'^, then uw G Wj"^'''. We estimate the product uw as 
follows. Let V G W^^^ and consider for any domain U with compact closure, vuw . By 
the above mentioned estimates, vu G W^''^{U). Thus the integral is well defined, and the 
inequality IJ^vuwl < C\\v\\w^,q(^u^\\u\\w2,pi^i/)\\w\\w-^.p{u) holds. But w G W^'''{U) 
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was arbitrary. It follows that uw defines a bounded linear functional on W^''^{U), and 
hence uw € W^^'^. 

Lemma 2.1 (Bianchi identities). Consider a Riemann manifold (M, hij) of dimension n, 
with metric hij G W^^f, p > n. Then the first Bianchi identity holds for Rijki- Further, 
the second Bianchi identity 

and the contracted second Bianchi identity 

V'^^, - iv,i? = 

are valid in the sense of distributions. 

Proof. The first two statements are clear from the product estimates and the definition of V 
and Rijki ■ It is most convenient to prove the Bianchi identity using the Cartan formalism. 
Let ^ be an orthomomal coframe. The structure equations are 

du) + CO Au) = Cl 

Assuming hij G W^^f we have cv € W^f and Q € L^. 

The Bianchi identity is the statement d^d'^d'^ = 0, where is the co variant exterior 
derivative. Recall that on a section of a tensor bundle, d'^d'^s = {du + ^[cj Acj])s = fls 
andona5o(n)-valuedtensor, such as f2, d'^H = dH + [to AH]. Evaluating d'^d'^d'^ gives 

rf'^Q = dQ + [ujAQ] 

= d'^w + ^d[LJ A w] + [w A n] 

Now, = on distributions. Further, expanding the other terms in the right hand side 
gives products of elements of W^^ and Lf^^. Therefore the standard algebraic identities 
hold to show that d'^'Q = in the sense of distributions. This is equivalent to the statement 
that ^[mRki]i = in the sense of distributions. Contracting this identity twice gives by 
the standard argument (making use of the product estimates stated above, and the fact that 
^mRkii^ € to justify the contraction), the identity 

V'Ri, - \v,R = 

which holds in the sense of distributions. □ 

For a domain f2 with boundary dO,, let xq denote the characteristic function of O, and 
trgo the trace to the boundary. The following Lemma characterizes the divergence free 
tensors supported on a domain. 

Lemma 2.2. Let (Af, h j j) he a Riemann manifold of dimension n with metric of class 
^\oc' P ^ ^- ^ bounded domain compactly contained in M. Assume that il has 

boundary dfl, and let Tij be a symmetric tensor of class W^^. Then 

^\TijXn) = (V'Ty)xn e 
if and only if the zero traction condition 

{tranTij)n^ = 
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holds in the sense of distributions, if and only if (V*ry )xsi = and the zero traction 
condition holds. 

Proof. Let e Cg^. Then we have 



where is the outward normal of dVl. This implies that with 1/^+1/^= 1, the inequality 



3. The field equations for a static, self-gravitating elastic body 

We first consider a variational formulation of a self-gravitating elastic body in a 3h-1 
dimensional spacetime (A^, and then specialize this to the static case. The body B is 
a 3-manifold, possibly with boundary. We use coordinates on B, and x°' on spacetime. 
In the Eulerian formulation of elasticity, the body is described by configurations f : Ai ^ 
B. The total Lagrangian density for the Einstein-matter system under consideration is, 
setting the speed of Ught c = 1 for convenience. 



where p = p{f, df, g) is the energy density of the materical in its own rest frame. General 
covariance implies that p is of the form p — p{f^,"f^^) where ~ f^,af^,i3g°'^^ 
where — daf^. Geometrically, this means the following: there is a function p on the 
bundle of contravariant, symmetric two-tensors over B, and p{f^, J^^) is the composition 
po where g"^ is the inverse metric and /, is pushforward under the map /, acting 

on contravariant two-tensors. For an equivalent description of elastic materials see lfT2l 
[13 1 . These references require in addition that B be furnished with a crystalline structure, 
which is essentially a choice of three linearly independent vector fields on B describing 
(the continuum limit of) the crystal lattice. More precisely, a crystalline structure is a linear 
subspace of the space of vector fields on B with the following property: for all points of B 
the evaluation map, which sends vector fields on B to tangent vectors at this point of B, is 
an isomorphism, when restricted to this subspace. Our assumptions, in section l33] below. 
will render such a choice of vector fields which Lie commute. This means there are no no 
dislocations in the crystal lattice. 

We now specialize to the static case. Let = K x M, where the space manifold M is 
diffeomorphic to Eucidean 3-space, M = M|. Further, we take the body S to be a bounded 
open domain in Eucliden 3-space, B C Mg. We refer to Mg as the extended body, and will 
use coordinate X"^ also on M^. The body B will be assumed to have smooth boundary 





holds if and only if {\iQ^T'^^)n^ — 0. This proves the Lemma. 



□ 
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dB, and the closure B U dB will be denoted by B. Letting (a;") ~ (t, x^) where are 
coordinates on M, we can write the static spacetime metric in the form 

(3.1) gapdx^dx^ = -e^^ dt^ + e-"^^^ h.jdx^dx' 

where hij depend only on x*. Further, the configurations f : M B aie assumed to 
depend only on x*. Assuming that a volume form Vabc, on B is given, we may introduce 
the particle number density n by 

(3.2) f^4x)f''^,{x)f,k{x)VABcif{x)) = nix) e,,k{x) 

where €ijk is the volume element of hij. Note that the actual number density defined with 
respect to the metric e^^^hij is e^^n. We assume the configurations / are orientation 
preserving in the sense that n is positive on f~^{B). 
Let 

(3.3) H^^ ^ f^^.f^^jh'^ 

„AB _ ^2U ttAB 



SO that = e^^ H^" . Note that equation (I3.2l i impHes 

Qu'^H^^'h^^'h^^'VabcVa'B'C' 
The Lagrangian density C is in terms of these variables, modulo a total divergence, 

(3.4) C^--^y/h{R-2\VU\'^)^e^ntVh 

Here, R is the scalar curvature of hij, |VC/p — h^^V iUV jU , and the relativistic stored 
energy function e, defined by p = ne, is of the form e = e(/'^, e^^ H^^), where e is a 
smooth function of its arguments. In particular, by the chain rule, we have de/dH^^ = 

Suppose that a non-relativistic stored energy function w(/'4,fs:'4'^) is given, for exam- 
ple one suggested by experiment, where K^^ is the non-relativistic analogue of H^^ . 
A relativistic stored energy function e corresponding to w can be defined as a sum of the 
specific rest mass e and the relativistic analogue w {f'^,e'^"H^^) of the stored energy 
function. The specific rest mass is defined such that IVabc is the rest mass distribution 
of the material in its natural state. It can be shown that if the dependence on the light 
speed c taken properly into account, the field equations tend to those of the corresponding 
Newtonian model when we let c /* cx). See |4| for details. 

3. L Field equations in Eulerian form. In order to write the field equations, we introduce 
the stress tensor a. We will need the form of the stress tensor on the body and in space. 
These are given by 

(3.5) (TAB = ~'^ gjjAB ' f^ij = nf^Af^ J (TAB, CFi^ ^ ^i asc 



We remark that our convention for the stress here is that used in standard nonrelativistic 
elasticity, as opposed to the usual one in general relativity. It is important to note that the 
elastic quantities such as H^^, Cij, viewed as functions on space, are only defined on 
f^^{B). The Euler-Lagrange equations resulting from the Lagrangian ( 13. 4t are 

(3.6a) Vj(e^a,J) = e^(ne - a,')V,C/ in f-\B), - 

(3.6b) A,,U = 4T:Ge^ine^c7i')xf-^s) in K| 

(3.6c) - 8^G(e,y - e^(T„- ^5 
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where 

(3.7) = -i^[v,C/V,C/- i/i,,|V[/p]. 

The system ( 13.61 ) is equivalent to the 4-dimensional Einstein equations 

(3.8) = SirCTf^iy, 

where Gp,y is the Einstein tensor of the static Lorentz metric given by (13.1b on R x M, and 

(3.9) Tf^^^dzf^dx" = e^ie^'^ne (u^rfa;^)2 - cr^^ dx'dx^ , u^d,, = e'^ dt 

The equations ( I3.6bl3.6cl ) together imply the elasticity equation ( I3.6al ). The reason is 
that the contracted Bianchi identity for Gy , V'Gy = 0, implies that the right hand side 
of (I3.6cl l has vanishing divergence, and in particular the divergence of the compactly sup- 
ported term aijXf-^(B) must be well defined. By Lemma l2!2l this implies that the zero 
traction boundary condition in equation ( I3.6al ) holds, and hence by equations ( I3.6bl ) and 
( 13.7b . equation ( I3.6ab follows. 

Let (5 be a fixed background metric on M, which we will take to be Euclidean, and let 
r* j. be the Christoffel symbol of 5. Then, with 

(3.10) F* = /i^'=(r},-fj,), 

— is the tension field of the identity map (Af, hij) {M, S), and we have the identity 

(3.11) i?„ = -^Ahh,, + V(,y,) + Qy (/i, dh), 

where Ahhij is the scalar Laplacian defined with respect to hij, acting on the components 
of hij and Qij is quadratic in dh. We use the notation ij^j) = |(ty + tji) for the sym- 
metrization of a tensor. In particular, hij i-^ Rij — V(iVj) is a quasihnear elliptic operator, 
while hij I— > Rij fails to be elliptic. This failure is essentially due to the covariance of Rij . 
It follows that also the system ( 13.61 ) fails to be elliptic. In order to construct solutions to 
( 13.6b . we will replace equation (I3.6cl ) by the reduced system which results from replacing 
Rij by Rij — V(jV^). The modified system which we will consider is of the form 

(3.12) - iA^% + Q,,{h, dh) = 2V,UV,U - 8^Ge^(ay - /ly 

3.2. Field equations in material form. In the Eulerian formulation above, the elastic- 
ity equation ( I3.6al ) is a nonlinear system with Neumann type boundary conditions on the 
domain f^^{B) which depends on the unknown configuration /. We will avoid dealing 
directly with this "free boundary" aspect of the system (13.6b by passing to the material, or 
Lagrangian form of the system. In this picture the configurations f : M ^ B are replaced 
by deformations, i.e maps (p : B ^ M satisfying (/> = f^^ on B. Recall that the body B is 
a bounded open domain in Kg, the extended body, and that B is assumed to have a smooth 
boundary dB. We will assume throughout the rest of this paper that B is connected. See 
remark l4!4l for discussion on this point. 

We assume given a diffeomorphism i : Rg K|, which is an isometry, i*S — (5^. In 
Cartesian coordinate systems X'^ on Mg and on R.^ where Sb and S have components 
Sab and Sij respectively, i can be assumed to have the form = 5\X^, so that 

iSai' = S\. Since B has smooth boundary, functions on B can be extended to the whole 
space, and in particular, given (p, there is an extension (p : M| M|, depending smoothly 
on (f>, such that (^{X) — i{X) for X outside some large ball. 
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In the material picture, the dependent variables /, U, hij are replaced by the fields 
0, U , hij which we now introduce. As mentioned above, (f) is assumed to be a diffeo- 
morphism B <t>{B) C M^, and the extension (f) of which depends real analytically on 
0, is used to define the fields U = U o (p, a function on Mg, and hij — hij o cj), a metric on 
Mg. We will use the same symbols for these fields restricted to B. 

Remark 3.1. It is important to note that hij ^ 4)* hij, since we are only pulling back the 
components of hij in the coordinate system (x*), not the tensor itself. In particular, hij 
does not transform as a tensor and is more propertly viewed as a collection of scalars. See 
appendixlAlfor discussion. 

The equation ( I3.2l i defining n can be written in the form f*V ~ nfih- where is the 
volume element of h. Defining J = n^^, we have (f)*iih — JV. The Piola transform of 
(Tj^ can now be written in the form 

With this notation we have in particular the relation V^o'i'^ = Ji\I jffi^) o (p. To derive 
the material version of ( |3.6at one may use this relation directly, or proceed by first pulling 
back the matter Lagrangian to B and then applying the variational principle, see IS]. One 
finds 

(3.13) VA{e^a,^)=e^[e-^]m in a/nA|as = 

n 

Here \/ A{e^ (Ji^) is defined in terms of the volume element V and does not involve a 
choice of metric on B. We have 

The bars in Eg. ( 13.131 1 correspond to the convention that f^^i be replaced by ip^i defined 
as a functional of (j> by 

(3.14) = S^B 
and H^^ be changed into 

(3.15) i?^^ = V^iV^j/^i, 

thus H^^ is the inverse of 4>*hij. Note ip"^, is defined on M|. With e = i{X, e^^ H^^) 
and H^^ understood in this sense, we have the identity 

(3.16) 



U<p\A 

In particular, di"^ — ip^ iGBcH'--^^, and hence ( 13.131 1 is a second order equation for (p. For 
barred quantities the corresponding rule for bars of derivatives gives for U 

(3.17) W^i'^^^AU 

For equations ( |3.6bl ) and ( 13.121 ) we simply replace each term by its barred version, i.e. 

(3.18) 'KhU ^ AT:Ge^ {nl - ai^)xB inM| 



Further, the covariance of the Laplacian implies A/jC/ — {/S.hU) o <p — Aj,^^(C/ o 0), so 
that in this expression, the puUback (f)*h appears. This tensor is given by the inverse of 
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H"^^. It follows that AhU involves second derivatives of (p. An analogous remark applies 
to 

(3.19) - ^Aj^ + Q,,(h,dh) = 2(V,C/)(Vjf/) - 8nGe^ia,,-h~ai')xB- 

3.3. Constitutive conditions. We shall need to assume that the elastic material is able to 
relax in Euclidean space (which in particular implies the absence of dislocations, see the 
remark at the beginning of section |3]l. Further, the stored energy function e must be such 
that the linearized elasticity operator is elliptic. In fact, we shall assume the pointwise 
stability condition. The constitutive conditions for e are thus formulated as follows. There 
should exist a Euclidean metric — Sab dX^dX^ on B (we will use the same symbol 
for this metric extended to Mg) such that 

(3.20) €{X) = e\^u=aM=5B) > C, ^ q^ab ^ \(u=o,H=Se) =0 in S 
and 

(3.21) Labcd TV^^iV^^ > C'{5ca5bd + 5cb5ad) iV^^^V^^ in B, 
where 

(3.22) LabCd{X) ■= ygffABQffCD J \{U=0,H=Sb) ■ 

and C, C are positive constants. The condition ( 13.21b is just the pointwise stability condi- 
tion ( 12.1b discussed in section IZ2] The quantity I appearing in ( 13.20b is the rest mass term 
in the relativistic stored energy function, as discussed above. For physical reasons, and 
in fact for hyperbolicity of the time dependent theory, it is necessary to assume that C is 
positive. However, for the purposes of the present paper, this condition could be dropped. 
We shall assume, for simplicity, that Vabc is the volume form associated with Sb (i.e. that 
^123 = 1 in Euclidean coordinates), so that n^/h — det{df ). 



4. Analytical setting 

We will use the implicit function theorem to construct, for small values of Newton's 
constant G, static self-gravitating elastic bodies near the reference state described in section 
13.31 We will use the field equations in the material form given by ( 13.13b . ( 13.181 ) and ( 13.191 ). 

We will now introduce the analytical setting where this work will be carried out. Fix a 
weight S E (—1,1/2). This choice of S determines the weighted Sobolev spaces which will 
be used in the implicit function theorem argument below. The range of weights for which 
the isomorphism property for A holds is (—1,0) but we shall need (5 e (—1, 1/2) later on. 
Further, we fix p > 3 to be used in setting up the function spaces which will appear in our 
argument. The body ;8 is a bounded open domain in Rg, the extended body, with smooth 
boundary dB. Under these conditions, the trace to the boundary tr^g is a continuous linear 
map W^'P{B) B^-^/P'P{dB), and there is a bounded, linear extension operator E : 
W^'P{B) — > W^^f{M.^), see the discussion in section ITTI or (T]. The spaces which will be 
used in the implicit function theorem argument are Bi = W'^'P{B) x Wg'^ x Eg'^, and let 
B2 = [LP{B) X B^-^/P'P{dB)] X L^_2 x ^5-2- Thus, Bi is a Banach manifold, and B2 
is a Banach space. 

The residuals of equations ( 13.13b . ( 13.181 ) and ( 13.191 ) define a map T : R x Bi ^ 
B2, T — T[G^ Z), where we use Z — {(p, U, hij) to denote a general element of Bi. 
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We assume that is a diffeomorphism onto its image. Thus, T has components J- = 
{Tff,, Tij , Th), corresponding to the components of B2, given by 

(4.1a) (^VA(e%-^)-e^[e--^]^, traal^r/w) 

(4.1b) Tv ~^T^Ge^{nl-al)xB 

(4.1c) Th = -^Ahh.j +Q,j(hM) - 2WfW^ + SttGc" {a,j -h~a/)xB 

Recall from the discussion in section [J!2] that the extension is needed in the definition of 
!F. The proof of the following Lemma is a straightforward construction involving the use 
of the extension operator E and a cutoff function, and is left to the reader. 

Lemma 4.1. Fix some Xq G B. There are constants /i > 0, i? such that for each (j) : B ^ 
^S' 1 1*?^" i| I VF^'^ — ^' ^here is an extension (j) : — > M.'^, which depends real analytically 
on (j). The extension <f> can be chosen such that the map i— > is given by a continuous 
linear map from W^'P{B) to W^'P{Br{Xo)), 4){X) = \{X)forX G M| \ Br{Xo), and 

r 1 e <f (K|). 

The equation to be solved is J-{G, Z) = 0. The material form of the reference state is 
given by 

Zo = (i,0,(5y oi) e Bi. 

The map JF defined by (14.1b is easily verified to satisfy J-{0, Zq) ^ and to map Bi — > B2 
locally near the reference state Zq. 

4. 1 . Differentiability of J-. In order to apply the implicit function theorem, we must ver- 
ify that T is as a map R x _Bi — > B2, in the arguments (G, Z), Z = (0, U, hij), near 
(0, Zo). In fact, is real analytic if the stored energy function e is real analytic) near Zq. 
It is clear the dependence on G is smooth. We will freely make use of the standard fact 
that if / is a smooth function, and u S W^'P{B), p > 3, then u 1-^ f{u) is a smooth 
mapping W^'P{B) W^'P{B), as well as the corresponding statement which holds for 
weighted Sobolev spaces. The map u f{u) is sometimes called a Nemytskii operator. 
We consider the dependence on Z for each term separately. 

It is straightforward to see that tp^i depends smoothly on (j). This means that in view of 
( 13.16b and the smoothness of the stored energy function, for JF^, we note that depends 
smoothly on Z. Expanding the definition of it is clear this depends smoothly on Z. 

Next consider J^jj. Note that AhU is the covariant Laplacian in the metric Hab, where 
Hab is the inverse of H^^ given by ( 13.151 ). and thus AhU depends smoothly on Z. Using 
the fact that — hai^cj)^ _a7 the lower order terms in JFy are seen to be smooth in Z 

Finally, we consider T^. The discussion above of AhU applies also to Ahhij. The 
quadratic term Qij{h,dh) is evaluated by replacing each occurrence of hij by hij and 
dkhij by 

In view of the regularity assumptions we find that Qij {h, dh) is smooth in Z. It is straigh- 
forward to analyze the remaining terms in JF/i in the same manner. We have now proved 

Lemma 4.2. The map T : Bi ^ B2 is near (0, Zq). 
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4.2. The Frechet derivative D2T{Q, Zq). Next we calculate the Frechet derivative 1)2^(0, Zq) 
and consider its properties. As we shall see, this derivative is not an isomorphism Tzg Bi 
B2, and therefore it will be necessary to modify the system of equations by applying a pro- 
jection before applying the implicit function theorem. 

We will denote by C, v, k the infinitesimal variations of the fields </>, U, hij. Using the 
notation of section [331 let 

The tensor S/!j ^ corresponds to the tensor Ai^ k of section l272l A calculation shows 

D^a,^{0,Z^).C^S,YdBC' 

D-^a,^{Q, Zo).k = 5//(5'"5Sfc™„ 

Let be a Killing field in the reference metric 6ij, in Euclidean coordinates = 
a' + P'jX', with a\f3'j constants, (3'^ = and define f (X) = CiK^))- Then we 

have SaC = f^^k^A' ^^'^ hence due to the antisymmetry of /3*fe, 

(4.2) ^Af^z^'' = 

Let now Sa^"^ denote any combination of the Frechet derivatives of a^"^, evaluated at 
(0, Zq). Assuming we use a coordinate system where 1^123 = 1, we have due to 
(I4.2I 1. and Stokes' theorem, the important relation 



= / CdAiSa,^) - / CiS^.^hA, 
Jb job 



OB 

where n"^ is the outward normal. This can be interpreted as saying that due to the natural 
boundary conditions, the linearized elasticity operator is automatically equilibrated at the 
reference configuration (0, Zq). It follows from the constitutive conditions stated in section 
13.31 that C I"* Dtp!Ftf,{Q, Zo)C, is elliptic, cf. the discussion in section IZ2] Therefore, the 
operator 

L>0j^0(O,Zo) : W^'P{B) ^ [LPiB) x B^-^^P^P{dB)] 

is Fredholm with kernel consisting of the Killing fields on B, and cokernel (in the sense of 
the natural pairing) consisting of the Killing fields on i(S), as above. 

The only nonzero contribution from the Frechet derivative of the first order term 

n 

in J^if, is 

J 

.V — —IdiV 

{0,Zo) 

This is a lower order term which cannot affect the Fredholm property of D2J-, but it should 
be noted that it is not a priori equilibrated, and therefore in general does not take values in 
the range of D^T^. By the above discussion we have 

i?2^0(O, Zo) = {dA{5a^^) - ed,v, traB{S&^^)nA) 

Considering the other components of D2J-, the only nonzero terms are the diagonal entries 
DfjTu and Dj^Th. These are given by 

DjjTu-v — Aw 
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and 

where A — S^^dAds is the Laplacian in the EucHdean background metric on M|. The 
operator A is an isomorphism ^5-2 ^ ^ 0)' l^l- 

It follows from the above discussion that the Frechet derivative D2J-{0, Zq) can be 
represented as the matrix of operators 

DrpTrj, DuT^ DhTfjA 
A 

-ia; 

(where the entries are evaluated at (0, Zq))- In particular, the matrix is upper triangular, 
and the diagonal entries are isomorphisms, with the exception for D^!F^{0, Zq) which is 
Fredhokn with nontrivial kernel and cokernel, cf. the discussion above. The off diagonal 
terms are bounded operators. Therefore, if we compose !F with a projection which in the 
first component maps onto the range of D^!F^{0, Zo), and restrict the domain of definition 
of to a subspace transverse to its kernel, the resulting map will have Frechet derivative 
at (0, Zq) which is an isomorphism, which will allow us to apply the implicit function 
theorem. The projection which will be used is introduced in the next section. 

4.3. Projection. Introduce the projection operator : B2 i?2, which acts as the iden- 
tity in the second and third components of B2 and is defined in the first component of B2 as 
the unique projection along the cokernel of D^!F^{0, Zq) onto the range of D^T^{0, Zq), 
which leaves the boundary data in the first component of B2 unchanged. We use the the 
label B to indicate that Pg operates on fields on the body and the extended body. We shall 
later need to transport the projection operator to fields on M|. 

In order to give the explicit form of the action of Pg in the first component of B2, let 
(bi, Ti) denote pairs of elements in LP{B) x B^~^/P'P{dB). By a slight abuse of notation, 
denote this operator too by Pg. We require that PB(6i, Ti) = {b'i, n), satisfying 

(4.3) / Cb\ - / en 

JB J as 

for all Killing fields Pairs Ti) satisfying this condition will be called equihbrated. 
Since Pg is a projection along the cokernel, b' i must be of the form 

b'i =bi- rjiXB 

where r/i = ai + PijX^ , with ai, /3ij constants satisfying /3y — —/3ji. We remark that the 
bilinear pairing (^,77) 1-^ ^'?7i is nondegenerate, it is simply the (B) inner product on 
the space of Killing fields. Further, the map 

C / Cb.. - / en 

JB JOB 

defines a linear functional on the space of killing fields, for each pair {bi,Ti). Therefore, 
there is a unique rji of the form given above such that b'i — bi~ rjiXB satisfies ( 14. 3t for all 
Killing fields f . 
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4.4. Existence of solutions to the projected system. We are now in a position to apply 
the implicit function theorem to prove 

Proposition 4.3. Let J- : Bi B2 be map defined by ( 14.71 ) and let Pg be defined as in 
section \4~3\ Then, for sufficiently small values of Newton's constant G, there is a solution 
Z ~ Z{G), where Z = (0, U , hij), to the reduced, projected equation for self-gravitating 
elastostatics given by 

(4.4) Pe^(G, Z) = 0. 

In particular, for any e > 0, there is a G > 0, such that Z{G) satisfies the inequality 

(4.5) 110- i||vi/2,P(z3) + 11% - '^uIIm/2,p + ||C7||^y2,p < e. 

Proof. Let Y denote the range of and let X be the subspace of Tz„Bi, such that {(f) — 
iy{Xa) = and 6'^ ^^c{A^B]{(t> - i)'(^o) = holds at some point Xq e B. We have 
already shown that T : x Bi B2 is and it follows from the definition of Pg that 
FbT : R X X Y is Since T{0, Zq) = 0, the Frechet derivative of Fs^iG, Z) 
with respect to Z, evaluated at (0, Zq) is Pg£'2.?^(0, Zq) which we will denote by A. It is 
clear from the discussion above that A is Fredholm with trivial cokernel. Therefore all that 
remains to be checked is that ker A is trivial. To see this, note that kei D^T^ consists of 
Killing fields, by the discussion in section |Z2l Prescribing the value and antisymmetrized 
derivative of a Killing field at one point determines the Killing field in all of B. This implies 
that A : X ^ Y has trivial kernel, and it is therefore an isomorphism. Thus, the implicit 
function theorem for Banach spaces |20| applies to prove the existence of solutions to the 
equation PbJF(G, Z) — 0, for small values of G. Since Z{G) depends continuously on G 
in B2, the inequality (I4.5l l follows. □ 

Remark 4.4. Up to this stage, the connectedness of the body B has not played at role in 
our arguments. In particular, in the case of a body with N > 1 connected components, 
there is a collection of N elastic fields cj) and an extension common to all of them. The 
equilibration argument of section |5] below fails, however, for the simple reason that we 
have one equation, namely ( 15.151 ). but N elastic fields. 

In the next section, we study the solutions provided by Proposition |43] and prove that 
they in fact represent solutions of the full system ( 13.61 ). 

5. Equilibration 

Let {(j), U, hij) e Si be a solution to the reduced, projected system ( 14.4b as in Propo- 
sition 14.31 and let be the extension of (f) provided by Lemma 14.11 Define U, hij by 
U = U o and hij — hij o (f)^^. Then U and hij are elements of VF^'^(K|) and 
Eg'^{M.^), respectively. Further, we set / = (p^^- Then / G W^f and / = outside a 
large ball. 

Let (j) be the extension of <j) provided by Lemma l4n and let U, hij be defined hy U = 
U o hij = hij o (j)^^. The following corollary to Proposition l4.3l is immediate in view 
of the composition properties of Sobolev functions. 

Corollary 5.1. For any e > 0, there is a G > such that the inequality 

(5.1) \\(l>-i\\w2.P(B) + - Sij\\yi,2,p + \\U\\^y2,p < e. 

holds. 
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5.1. Eulerian form of the projected system. Next we introduce the projection operator 
in space which corresponds to Pg. From this point on, we are only interested in the action 
of Pg in the first component of Bi. We will be dealing with solutions to the projected 
system PgJF = 0, in particular this means that the boundary condition (Ti"^ n^|^ ^ = will 
be satisfied, due to the nature of the projection operator defined in section |4~3T Therefore 
it is convenient to write Pg6 = Pz3(6, 0). 
Recall the change of variables formula 

u{X)dV{X)= i uo f(x)n{x)dfih{x) 
Jf-HB) 



which results from the definition of n in equation ( 13.21 ). This implies 



eiXMX)dV{X) ^ / e°.fix)n{x)b,of{x)dfih 
Jf-HB) 

Letting bi = "S/ a{g^ — e'^(e — ^) djU, we have, using the properties of the Piola 
transform, 

Based on this, we define the space version Py-i(g) of the projection Pg by 

P/-i(B)(n.(6o/)) = n(PB6)o/ 

Pg is aprojectionby construction, and it follows from the definition of Pj^-i(g) that it also 
is a projection. Due to this definition, the equation PgJF^ = Oimplies P^-i(g')(njr0o/) — 
0. This relation is written more explicitely as equation ( I5.2al i below. Summarizing, the 
triple ((/), [/, hij) we have constructed satisfies the following set of equations, 

(5.2a) 

^f-'{B) (Vj(e^a,^") - e^(ne - ai')d,U) =0 in r\B), 
(5.2b) =0 indr\B) 



(5.2c) AhU = A^Ge"{n€-ai')xf-i(B) 



inlt 



■s 



(5.2d) G„- - (V(,y,) - -h^,^lV') = 87rG(e,y - e^tr,, Xf-^(B)) in M| 



1 
2 

Recall that equation (I5.2db is equivalent to the space version of equation (|3.19l l. 



(5.3) - ^Aft/iy + Qij{h,dh) ^ 2V,[/Vj[/ - 87rGe^(CTy - h,jcri')xf-^B}- 

Since the body is bounded, the right hand sides of equation ( I5.2cl i and the contribution 
from the stress in equation ( 15.31 ) have compact support, and hence U, hij satisfy the reduced 
vacuum static Einstein equations near spatial infinity. This implies that we can get more 
information about their asymptotic behavior than is a priori given from the implicit function 
theorem argument. 

Lemma 5.2. U, hij are of the form 

(5.4) - ^ + Ui2) 

la 

(5.5) h^j = % + -p- + h(2) ij 
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for constants mu, "fij, with ^(2) , t^(2) G ^25^- Further, for sufficiently small G, we have 
the estimates 

(5-6) \\U(2)\\w^,p + \mu\ < C{\\h,j - S,j\\^2.„ + \\<p - i\\w^.v(B)) 

(5.7) ||/l(2) ^j\\w^^P + II7II < C{\\h,j - Sij\\^2.p + \\(j)-i\\w^.P(B)) 

where \\j\\ = (E«j7y-)^/^- 

Remark 5.3. If hij, U were solutions of the full static vacuum Einstein equations near 
infinity, with the same falloff conditions, one could conclude that the 0(l/r) term in hij 
vanishes, following the argument in [17J . see also [7J. However, in our situation U, hij 
satisfies the reduced static vacuum Einstein equations near infinity and this argument does 
not apply directly. 

Proof. By construction, we have hij — Sij G W^'^ and U G W^'^, with (5 e (—1,-1/2). 
This, together with equation (I5.2cl ) implies using the product estimates (12.9b to expand 
A,,C/, that 

with 

< C(||ft,ij - 5ij\\^y2,p + 110 - i||i4'2,p(f3)), 

where in the last step we estimated ne — cr;' in terms of hij, (j) — i. 

We follow the proof of ID Theorem 17], see also the discussion in section |23] Let u 
satisfy 

AeM / e LP_5_2 

By the isomorphism property of Ae : W^'^ ^?-2' '•^^ inclusion ^2(5-2 ^5-2' 
we have 

Fix some large radius R, and a point G K."^. Let Bji = Bj^{xo), and let Ej^ = R"^ — 
Bji. By assumption, d G (-1,-1/2), so that —2 < 2(5 < —1. In particular 26 is 
nonexceptional, and hence Ae : — > ^25-2 Fredholm. Hence there is a u G W^2i^ 
such that 

Ae(M — w) =0 in En 

It follows that u — V satisfies 

M- w = — + 0(l/r2) ini;^? 



for some constant m. Here the term 0(l/r^) comes from the expansion of harmonic 
functions, and is in W^f ■ Let C, be 
Bn{xQ) and with |aCI < GjR. Set 



functions, and is in W'^f ■ Let C be a cutoff function such that ^ = 1 in i?27?, C = 



m 

w — u 

r 

Then we have w = uin Bj^{xo). Further, we have 



Aeu; = /-A,(-C) 
r 

where Ae(— C) has compact support and we can estimate 

\\^^(jO\\Ll,_,<C\m\ 
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Since 

where F = 1/r is the fundamental solution of Ag we have 

H<c||/|U.(K3) <qi/||L5,_, 

Here we used the inclusion ^25-2 -^25-2 ^-3' ^^'^ ^'^^^ ^^^^ = ^-3 ™ ^ 
dimensions. Thus we have 

||AeHlLg_<qi/llL5_ 

Using the fact that ker Ag is trivial in W^f together with the Fredholm property, we now 
have 

IHIw-<c||/ILg_ 

This argument proves ( 15.41) and ( 15.61 ). Applying the same argument to hij shows 

with constant and /i(2)ij € ^2s' ^ ™d gives an estimate for h(^2)ijjlij in terms of 
hij — Sij, U, and cj) — i. This gives iS.li . since we have an estimate for U in terms of 
hij — Sij, ip — i from ( 15. 6t . □ 

5.2. The tension field. The next three lemmas state the properties of the tension field V^, 
defined by ( 13.101 ), which we shall need. We may assume r*^. — 0. Let the operator L be 
defined by 

where A/j, i?y are the covariant Laplacian and the Ricci tensor of hij, respectively. 

Lemma 5.4. Assume we are in a Cartesian coordinate system so that Sij — Sij and Tij — 
0. Let r — \x\ and set 9^ ~ jr. Let hij, jij be as in ( 15.51 ). Suppose LV G -^5_3. where 
V is the tension field ofh. Then ifG is sufficiently small, V € W^Jl'i i^nd the estimate 

(5.8) ll^llwl- <C||iV^||Lj_3 
holds. Further, V is of the form 

(5.9) ^ " 2^ ^ ^'"™^') + ^(3) 
where V^^^ G M^a^i^ p andj'-j = S''''jkj- 

Proof. The operator L is an elliptic system, which is of diagonal type. Therefore the results 
discussed in section |23] apply to L. In particular, L is asymptotic to A with rate S. Thus, 
L : Wgf-^ ^5-3 is Fredholm. It follows from the definition of V and the assumptions 
on h that V G W^^l^i. If LV G L^_^, then elliptic regularity gives V G M^ll^i. For small 
data L has trivial kernel in this range of spaces, since L is then a small perturbation of Ag 
which has trivial kernel in W^'^, r < 0. By Corollarv l5.1l for G sufficiently small, the 
required condition on L will hold. Therefore, if G is sufficiently small, an estimate of the 
form (15.8b holds. 

Recall from Lemma |5?2l that 

^ij 

hij = Sij + —— + h^2) ij 
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with /i(2)ij e W^f- Let 



(5.10) = " 2^ (^'j"^' ^ '^'"™^') ■ 

Under our assumptions, T^j = 0. From the definition of V'', a calculation shows 

in particular V'' — H'^ = o[l/r'^^^^). The operator L is asymptotic to Ag with rate 5, see 
section l23] for this notion. It follows by an argument along the lines of the proof of Lemma 
lUthat 

yk ^ffk^ yk^ 

with y^gj £ ^25-1- '^i^ completes the proof of the Lemma. □ 

Let C be a Euchdean KilHng field of the form f (X) = f (i(^)) with f (x) ^ a' + 
P'^jX^ for a% constants such that = —f3^i. 



Lemma 5.5. Lef V be as in ( I5.9I >. TTien f/ie identity 

(5.11) y f o = y" o 

Proof. First note that both sides of equation (15.1 11 1 are well defined. The integrand on the 
left hand side has compact support by equation ( 15.151 ). Recall that / equals the identity 
map outside of a compact set. We have 

and 

V, = 0(l/r2) 

Therefore the integrand in the right hand side is an element of LFg_^ C L^. In order to 
justify the partial integration, let Bj^ = {\x\ < R} and 5;^ = {|.t| = R}. Then ( 15.1 lb is 
equivalent to 

lim / Cofmdfih^ lim / L{C o f)V,dtih. 
Gauss law applied to the integral over Bji gives 

(5.12) f C°fLV,dtJih^ I L(^'of)Vdiih 

J Be -I Br 

+ f Cof^kV^n'^dAh- f Wk{eof)n''V,dAh 

JSr J Sr 



Let 9 — X jr. With the asymptotic behavior of hij from Lemma 15721 we have n = 
9^ + 0(l/r) and r*j, = 0(l/r^). Further, the area element induced on Sh from hij 
differs from the standard area element of 5^ by a term of order 0{\/R). Recall that by 
Lemma |54l V*" = + V^^^ with H'' given by dSlOl i and V^:-^^ G W'^f^^. Since we are 
interested in the integral over Sr in the limit R oo, this shows that the only important 
terms in the boundary integrals in ( I5.12l i are 

(5.13) f p'',x'diHk9'dAR~ f p\W9''dAR 
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We consider the first of these terms. Thus let Z ~ f3''ix''diHkO\ the integrand in the 
first boundary integral above. At this point, we will be doing all calculations in the back- 
ground metric 5ij, and therefore it is convenient to lower all indices using Sij and sum over 
repeated indices. A calculation shows 

dlHk = [27ifc - IrnmSlk - 3{2jjk0l6j - JmmOlOk)] 

This gives using OkOk — 1, 

1 

Note that the second term in the right hand side vanishes due to the fact that Pki — 
—f3ik, but not the first. To show this vanishes we proceed as follows. We have dAn = 
R^dA where dA is the area element of the unit sphere S, so we can write J^^ ZdAn = 
ZdA. Recall the identity 

0keidA = ^5ki . 



Z = (ikiXidiHkOi = -—f3ki0ii-4'^jk0] + '^Immdk) 



This gives 

/ ZdAu = —Pkii-'ijjkSij + 2jmrn6ik) = 

JSr 

due to the fact that (3ki = —Pik- The integrand in the second term in ( 15.131 1 is 

which can be handled using the same method. It follows that 

lim ( [ Co f^kV'n^dAh - I Vfc(r o f)n''VdAr}\ = 
^'^•^ \Jsr Jsr J 

which completes the proof of the Lemma. □ 
Lemma 5.6. For sufficiently small G, there is a constant C such that the inequality 

^ij + II'/' l|llV2.P(B) j I^IIIV'll^yl^P^ 

holds, where for C = + (3^jx\ we write |^| = |a| + 
Proof. Using Lemma |53] we consider 

L(f o f)V,dfih 



We need to estimate ||L(^' o /)Vi| |li. We expand out o /), dropping the reference to 
/ for brevity, and writing ★ to denote a general contraction 

Using the form of hij given in equation ( 15. St . we have 

r^O{l/r^) + v, withweL^^^^ 

dr = 0{l/r^) + z, withzeLP_5_2 

R = 0{l/r^) + w, withweLPs_2 

Further, using the estimates from Lemma |5.2| the coefficients in the 0(l/r^) and 0{l/r^) 
terms can be estimated in terms of |7|, where 7^ is the constant metric in ( 15. 5t , and the 
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lower order terms can be estimated in terms of /i(2) y ■ The term involving d^^ has compact 
support and can be estimated in terms of — i. Multiplying o / by V, using — 
0{r) and V G VF^v?l, we find the result gives terms in -L^ g and L^g_2 coming from the 
0(l/r^), 0{l/r^) terms and the Sobolev terms, respectively. Since L^g_2 -^5-3 
C L^, this means that using Lemma l572l we have an estimate of the form 

\mof)v\\Li < c (i7i + \\h(2)^^\w^^■^^ + \\<p--A\wi^) n \\y\\wt\ 

< C (^\\h,j -Sij\\^2.p + \ \(t)-i\\w2.p(^B)) 1^1 1 1^1 1 w,";^^ 
This completes the proof. □ 

We now take, in the sense of distributions, the divergence of equation ( I5.2dl ). i.e. the 
equation 

(5.14) G,, - (V(,y,) - ^h,,ViV') = 87rG(e,, - e^a,, Xf^HB)) , 

which holds in R|. For the first term on the left in equation ( 15.141 1 we use the Bianchi 
identity from Lemma IZTI to conclude V*G,;j = in the sense of distributions. In order to 
take the divergence of the second term on the left, note that it follows from the definition 
of V that V G W^^l^i • The regularity assumptions on hij imply that the identity 

holds in the sense of distributions. For the first term on the right of ( 15.14b . we use equation 
(I5.2cl l. In the second term on the right in ( 15.14b . using that the boundary condition (I5.2bb 
is satisfied, we apply Lemma lzSl to get the identity 

In view of the projected elastic equation ( I5.2ab . this is equivalent to 

- V^(e,y - <T,jXf-^B)) 

where ly-i(g) is the identity operator in the space L'P{f^^{B)). Thus, defining the operator 
Lby 

we finally conclude that the equation 

(5.15) m = 16^G(I/-i(B) - P/-i(B)) [Vj(e^a,^") ~ e^(ne - ai')d,U] Xf-^s) 
holds in R|. 

We observe that equation (15.15b implies that LV is supported in f^^{B). Further, in 
view of the fact that Pj-i(j3) is a projection, we have 

(5.16) ff-i(B)LV = Q. 

Recall from section |4]that Pg was defined in terms of the Killing fields — 

with = a* + jX^ for a*, (3^j constants such that = —(3^1. It follows from the 

definition of Pj-i(b-), that we have 

0=/ {C ° f{x))^f~i{B)ZiXf-HB)diJ-h 
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for any ^' of the above form and any zi G ^^^^.(Ml). Now define the hnear mapping 

Q : Lg_3(M|) M^ by setting 



where {C(k)}k=i forms a basis for the space of Killing fields. We have the estimate 
Lemma 5.7. There is a constant, depending on f and hij such that the inequality 

(5.17) ||z||ip(y-i(B)) < C{\\Ff-^B)z\\Li-(f-^B)) + IIQ^X/-i(b)IIr'5) 

holds. 

Proof. From the definition of Py - 1 , we have 

2X/-i(i3) = P/-i(e)2X/-i(B) + "-(C ° .f)Xf-HB) 
where ( is a Killing field, and from the definition of Q, we have 

Q(k)ZX/-i(B) = / ('?(k) ° D^iC," f)Xi-HB)dl^h 



JB 



for a basis of Killing fields {S,(k)]^^=i- In view of the non-degeneracy of the pairing 
on the space of Killing fields on B, this means that there is a constant C such that | \n{C, o 
/)l lLp(/-i(i3)) ^ C*! |Q"-(C ° /)X/-i(8) I Ik". The inequalitv (15.17b now follows after an 
application of the triangle inequality. □ 

Proposition 5.8. For sufficiently small G, there is a constant C such that the estimate 
(5.18) \\V\\^.._,^<C\\QLV\W. 

holds. 

Proof. Due to the boundedness of B and f^^{B), there is a constant C such that 

lkX/-i(f5)llL?_3 < C'||z||LP(/-i(e)) 

Since V solves ( 15.151 1. LV is supported in f^^{B) so the above inequality applies to 
LV. Now making use of (15.161 1. the estimate (15.17b . together with the estimate (15.8b from 
Lemma l54l we find that for hij sufficiently close to Sij, i.e. for G suffiently close to 0, the 
estimate ( 15. 18b holds. □ 

5.3. The main theorem. We are now able to prove 

Theorem 5.9. For sufficiently small values of G, the solution to the reduced, projected, 
system of equations for a static, elastic, self-gravitating body, equations ( 15.21 ) is a solution 
to the full system ( li.6b of equations for a static, elastic, self- gravitating body. 

Proof. Recall the definition of Q. We have 

For small data, i.e. for sufficiently small G, ^ is an approximate Killing field for hij. 
Therefore by Lemma |531 there is for sufficiently small G a constant G such that 
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By ( 15.18b . this implies that after possibly decreasing G, there is a constant C such that 

< C' (\\h,j - (5y||^2,p + ||(/'- i||M/2,P(e)) llV'llw^.p 

0—1 \ / 0—1 

Hence, by Corollarv l5.ll for sufficiently small values of G, we have 

1, 



/2,P 



< 



Thus, in fact V = 0. Hence equation ( I5.2dl ) implies the full Einstein equation ( I3.6cl) . 
Further, by ( 15.15b we have that 

- P/-i(s)) [Vj(e^a,^') - e^ine - ai')d,U] Xf-HB) = 
while by equation ( I5.2al ). 

P/-i(6) (Vj(e^(T,^) - e^(ne - ai')d,U) = in r\B) 

It follows that equation ( I3.6ab is satisfied, and hence that the solution (</), U, hij) to the 
projected, reduced system, which was constructed using the implicit function theorem, is 
a solution to the full system ( 13.61 ). □ 



Appendix A. 

The way the spatial metric hij is treated in section |4] may seem surprising at first. We 
start by reducing the Einstein equation ( I3.6cb using harmonic gauge, and then pull back the 
components of hij in that gauge. Geometrically it might seem much more natural to start 
by pulling back the metric itself under Doing this would result in replacing ( I3.6cl ) by 

Gab = SirGiQAB - e^^aABX^) in 
where Gab is the Einstein tensor of Hab- Furthermore, we have that 

Qab = -^[DaUDbU - IhabH^'^'DcUDdU] 

and aAB is the extension of <tab- But then the quantity ai^ in equation ( 13.161 ) depends on 
only through ^p^i, since the remaining factors in 

(A.l) af = ^''^JaBcH^'' 

depend only on Hab- Let us now look at the nature of Eq. (13.13b . when (0. U . Hab) are 
used as the basic variables rather than {(f), U, hij). There are potentially two terms depend- 
ing on second derivatives of 0. One comes from taking a derivative of ip^i in equation 
dA.lb above. The other comes from writing the Christoffel symbols of hij entering the left 
hand side of (13.131 ) in terms of (f>. Remarkably, both these terms cancel! For the lineariza- 
tion of the residual of equation ( 13.131 ) at {4> — i, U ~ 0, Hab = ^ab) the only surviving 
term will involve first partial derivatives of Hab, i-e., terms which break diffeomorphism 
invariance on Rg. However note that the very presence of the a priori given domain B 
(which of course is the reason why this whole manoeuvre was performed) breaks diffeo- 
morphism invariance. It is thus no wonder, that, after going over to the material picture in 
the sense of using the full pull-back of hij, the elastic equation gives us terms which break 
diffeomorphism invariance. 
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Appendix B. 

The following Newtonian version of the argument in this paper was contributed by the 
anonymous referee. In the Newtonian theory the equations, with the standard definition of 
the stress, taking the specific mass to be 1, are 

(B.la) Vj(7i^ + nd,U = in .r\B), 

(B.lb) cri%j|g^_i^gj = 0, 

(B.lc) AC/ = 47rGnx/-i(g) in R| 

The argument to show that the projected version of dB.lal l impUes dB.lal) would go as 
follows. We define Xi to be the solution of 

(B.2) AX, = Vjcr,^' + nd,U 

tending to zero at infinity. Since 

Pi = Vjcr.j-' + ndiU 

the right hand side of ( IB.2b . has compact support, Xi has a multipole expansion in a neigh- 
borhood of infinity. 

Mi xWi. , 

Here 

Dij = / x-' pi{x)(fix. 
Jr3 



Writing 



where 8^^ ~ 8^ is given by (13.7b with hij replaced by dij, using the fact that Qij 
0{r^^) at infinity, and taking into account (IB.lbb . we obtain 

Af = 

= - / (a,^ + Qe)d?x 



Thus, 



where = Dji, in agreement with (15. 9b . Lemma 153] then applies with A in the role 
of L and Lemma |53] applies with 5ij in the role of hij. Proposition 15 . 8 1 also applies. The 
argument of the proof of Theorem l5.9l then applies to show that Xi ~ 0, therefore equation 
(iBlab holds. 



X, = -^ + 0{r-'), 
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